Small-amplitude normal modes of a vortex in a trapped Bose-Einstein condensate 
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We consider a cylindrically symmetric trap containing a small Bose-Einstein condensate with a singly 
quantized vortex on the axis of symmetry. A time-dependent variational Lagrangian analysis yields 
the small-amplitude dynamics of the vortex and the condensate, directly determining the equations 
of motion of the coupled normal modes. As found previously from the Bogoliubov equations, there 
are two rigid dipole modes and one anomalous mode with a negative frequency when seen in the 
laboratory frame. 
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I. INTRODUCTION 

After the achievement of Bose-Einstein condensates in 
trapped atomic gases , interest in these systems has 
increased dramatically. An especially intriguing ques- 
tion concerns condensates containing a vortex. The low- 
energy normal modes are important in characterizing the 
behavior of such a system. They have been studied exten- 
sively f|-[ll| , mostly with the Bogoliubov equations ]T^ ] 
that involve quantum-mechanical normal-mode ampli- 
tudes. Consequently, the associated dynamical motion 
of the vortex core and the condensate center of mass can 
only be obtained indirectly. Here, we use a variational 
Lagrangian procedure to provide a more direct and intu- 
itive treatment of the small-amplitude normal modes of a 
small (weakly interacting) Bose-Einstein condensate con- 
taining a singly quantized vortex. As in previous applica- 
tions |§,[f|, this approach yields a clear physical picture 
of the dynamics. We show that it reproduces the three 
lowest modes as found from the more intricate Bogoli- 
ubov analysis, thereby clarifying the physical interpreta- 
tion. This description is particularly interesting for the 
anomalous mode that causes the local instability of the 
vortex in a trap at rest fl5p|]i"o|] . 

The article is organized as follows: in Sec. II we derive 
the effective Lagrangian for a slightly off-center vortex 
in a small condensate. In the next section we solve the 
corresponding equations of motion for the vortex and the 
condensate and discuss the physical interpretation. Fi- 
nally we summarize our results in the conclusion. 



II. EFFECTIVE LAGRANGIAN 

We consider a small Bose-Einstein condensate with N 
particles containing a singly quantized vortex on the axis 
of symmetry, and observed in a frame rotating with an 
angular velocity Q. The harmonic trap is axisymmetric, 
with radial and axial frequencies loj_ and lo z (the ratio 
A = lo z /lo±_ characterizes the axial asymmetry). The s- 
wave scattering length a (assumed positive) characterizes 



the interparticle interactions. For a general condensate 
wave function \P, the Lagrangian for such a system is 
given by 



L(*) 



dV 



- (**9 t * - I'd***) 

<T (H - flL z ) * - 2vr7|*| 



(1) 



where Hq is the Hamiltonian for the noninteracting con- 
densate 



Ho = l 



-dl-dl + { X *+y*) + \(-dl + Z i) 



, (2) 



and L z = —i(xd y — yd x ) is the z-component of the 
angular-momentum operator. All quantities are ex- 
pressed in dimensionless units (the radial and axial co- 
ordinates are scaled with the radial and axial oscillator 
lengths d± = ^h/Muj± and d z = y/h/Mw z , and fre- 
quencies are scaled with the radial trap frequency u>±). 
The condensate wave function \& is normalized to unity 
and 7 = Na/d z is the small interaction parameter. 

Next we construct a trial condensate wave function. 
The condensate wave function is assumed to be un- 
changed with respect to the noninteracting case along 
the axis of symmetry, which allows us to use the ground- 
state gaussian <po(z) = n^ 1 / 4 exp (— z 2 /2). The ra- 
dial wave function is based on the vortex ground state, 
namely a two-dimensional gaussian times the vortex fac- 
tor x+iy = re %v . We are interested in the relative motion 
of the vortex which requires introducing time-dependent 
parameters for the vortex position. In the weak-coupling 
limit, the radius of the vortex core is comparable to 
the radius of the condensate, so that the displacement 
r o(t) = [xo(t),Vo(t)] of the vortex, the displacement 
rx(t) = [xi(t),yx(t)] of the condensate, and the induced 
velocity of the condensate a(t) — [a x (t),a y (t)] must all 
be included. Hence we use the following trial function 



ty v (x,y,z,t) 



x e 



C 

1 [(x-x 1 ) 2 + (y-y 1 ) 2 +z' 2 ] i( a:c x+a y y) 



[(x - x ) + i(y - y )] 



(3) 
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where C~ 2 = 1 + |rj — ro| 2 . Since we focus on the dynam- 
ical motion of the vortex and the condensate, we ignore 
the possibility of monopole or quadrupole modes, and 
the dimensionless gaussian widths are taken as one. The 
same trial function, Eq. (||), also serves to characterize 
the stability of the vortex in a small condensate Jl4| . 

In evaluating the integration in L(JS V ), we retain all 
terms up to second order in the time-dependent param- 
eters, yielding 
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This Lagrangian takes a simpler form when written in 
terms of new variables 8 = ri — ro and e = 2r\ — t*o: 
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contrast, the e mode with = represents a rigid dipole 
oscillation with t*o = r% (so that the vortex and the con- 
densate move together). 

The physically relevant solutions of the equation of mo- 
tion must represent quantum states with positive normal- 
ization. A convenient condition to check this normaliza- 
tion is 



/ dVi(n'*& -*'V) = -^, 



(11) 



where n' and $' are the complex fluctuations of the 
particle density and the velocity potential respectively. 
The density fluctuations are obtained by expanding 
the density n — \^ v \ 2 = no + Sn, where no — 
(r 2 /7r 3 / 2 ) exp(— r 2 ) exp(— z 2 ) is the unperturbed density 
for the condensate with a vortex. We find 

TV 



2tiq r 



(12) 



and use of the normal-mode amplitudes yields the ex- 
pression in plane polar coordinates 



Sn 



2n a { cos ip[r{e x - S x ) - r (e x - 25 x )] 

+ shap[r(e y - S y ) - r" 1 ^ - 2S y )]} . (13) 



Given the explicit form of the displacements S and e (see 
below) , this real expression can be converted to complex 
density fluctuations through 



Sn 



n e 



■n'*e iuJt . 



(14) 



The complex fluctuations in the velocity potential fol- 
low by expanding the phase of the condensate wave func- 
tion S — arctan[(y — yo)/(x — So)] + a ■ r ip + 5S, since 



S$ = {h/M)SS = &e 



(15) 



III. NORMAL MODES 

The effective Lagrangian leads to six coupled differ- 
ential equations for the displacements and the induced 
velocity. Elimination of the latter through 



e + z x (SI e — 8) 



(6) 



yields two uncoupled pairs of homogeneous equations for 
the displacements 8 and e: 



1 - 
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-sij s y (t) + s x (t) = o, 

-Si) S x (t) + S y (t) = 0, 



(SI 2 - l)e x (t) + 2Sle y (t) - e x (t) = , 
(SI 2 - l)e v (t) - 2Sle x (t) - e y (t) = . 



(7) 

(8) 

(9) 
(10) 



The 8 mode with e = will be seen to correspond to 
the anomalous mode; it has ro = 2t\ (so that the dis- 
placement of vortex is twice that of the condensate) . In 



The vortex part 5S V of this phase perturbation is SS V = 
(—yo x + x oy)/ r2 , and condensate- velocity part 5S a = a r 
can be rewritten with Eqs. @. In terms of the normal 
mode variables 8 and e, we find 



8S 



COS if 



sin ip 



r( , 



^ + rSl 
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and this expression can readily be converted into the 
complex velocity potential <&'. 

The equations of motion for the anomalous mode 8 
involve the characteristic frequency 



7 



SI. 



uj„ = — 1 

2V2^ 

and the solution has the form 

8(t) = -^=<5 (cosuj , - smu a t) , 



(17) 



(18) 
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where 5q is an infinitesimal amplitude and the factor 
l/y/2 fixes the normalization of the perturbation ampli- 
tudes (see below). Since e = for this mode, the real 
small perturbation in the density becomes 



5n a = —2no yr J (cos (p5 x + sin if 5 y ) , (19) 

and the complex density fluctuation (with time depen- 
dence e _i " at ) is 



So 



(20) 



As expected from the physical picture of a single-particle 
transition from the vortex condensate with unit angular 
momentum per particle to the unoccupied axisymmetric 
gaussian ground state, the quantum number of the z- 
component of the angular momentum in this mode is 
m = — 1. The real phase fluctuation for the anomalous 
mode is 



SS a = ( — h r ) (cos ipSy — smtp5 x ) 



with the corresponding complex velocity potential 



Sp ft 



(21) 



(22) 



For this anomalous mode, it is easy to verify that the 
normalization integral ( [TT| ) has the proper positive value. 
Since e = 0, Eq. ( |l8| ) shows that the amplitudes of the 
vortex r and the condensate r\ become 



r = V2 So {— cos uj a t , sin uj a t) , 
r\ = i\/2 <5 (- cosw a t , sinuj a t) 



(23) 
(24) 



they both execute small-amplitude in-phase motion with 
ro twice as large as r\. For = 0, the anomalous 
frequency (|l7]) is negative, and the motion is counter- 
clockwise in a positive direction (namely right-circular 
motion with positive helicity). With increasing rotation 
frequency f2, the negative anomalous frequency increases 
towards zero and vanishes at the metastable rotation fre- 
quency Sl m = 1 - 7/(2 v^7r) 0]. For Q > Q m , the 
apparent motion becomes clockwise (namely left-circular 
motion in a negative direction). The expressions found 
here for the density and velocity potential fluctuations n' a 
and <f>' a coincide with those found from the Bogoliubov 
approach || . If we had omitted either the induced veloc- 
ity or the displacement of the condensate, we would have 
found different frequencies, just as in the discussion of 
vortex stability jwjj. The present Lagrangian treatment 
has the advantage of providing a physical picture of the 
motion, clarifying the nature of the anomalous mode. In 
the weak-coupling limit, the vortex and condensate move 
in phase, in contrast to the behavior found in the strong- 
coupling (Thomas-Fermi) limit 



In a similar manner, we can treat the dipole modes, 
which have S = 0, so that ro = r% and the vortex moves 
rigidly with the condensate. Equations (||) and ([h]) show 
that the x and y motions are uncoupled if Q = 0. In gen- 
eral, however, there are two coupled normal modes with 
frequencies uj+ = 1 — Q and c<j_ = l + O. The notation w± 
reflects the general construction of the Bogoliubov am- 
plitudes u± = alj_\l/ and v± = a^^* for the dipole modes 
in the laboratory frame, where \t is any solution of the 
Gross-Pitaevskii equation in an axisymmetric trap; the 
resulting density perturbation is n'± oc e iv> e~ tu±i , see 
Eqs. (|28| ) and (|o|) below. It is easy to see that the gen- 
eral solution is a linear combination of the two distinct 
normal modes 

e + (t) = e_|_(0) (coscj+i , sinw+t) , (25) 
e_(t) = e_(0) (cosw_t,-sinw_i) , (26) 

where e_|_(0), e_(0) are infinitesimal small amplitudes. 

We can now check the normalization for the dipole 
modes. The first-order contribution to the density for 
both modes is given by 



Sn± = 2no ( r ) (cos <p e x ± + sin (p e y ±) , 



and hence the complex density fluctuations are 



n'± = e±(0)n [r e 



1 



(27) 



(28) 



For the phase fluctuations we find 
1 



SS± = 



(sin ip e x — cos tp e y ) , (29) 



and the resulting complex fluctuations for the velocity 
potential are 



± 2iM 



(30) 



independent of the external rotation frequency Q. With 
these expressions, the normalization integral is one for 
both normal modes with frequencies lo± . 

The final expression for the motion of the vortex and 
the condensate in the dipole modes \r(t) = ro(t) = ri(t)] 
is 



r±(t) = ro(coso;±f, ±sinw±t) 



(31) 



with an infinitesimal amplitude ro. We again find cir- 
cular trajectories, and here they are the same for the 
vortex and the condensate motion, confirming that the 
modes are indeed the rigid dipole oscillations of the vor- 
tex and the condensate. For stability reasons, only the 
region | \< 1 is relevant; thus the motion is always 
counter-clockwise (positive) for oj+ and clockwise (nega- 
tive) for uj- . Comparison with the Bogoliubov results || 



3 



shows exact agreement for O = 0. The rotation splits 
the dipole frequencies by twice the external rotation fre- 
quency. This result can be visualized by expressing the 
dipole modes (for O = 0) in terms of circular quanta, 
namely as circular motions in the positive and negative 
sense. For zero external rotation, the angular velocities 
are degenerate. In a frame rotating with O in the positive 
sense, the positive mode rotates at the difference 1 — 
and hence slows down. The negative mode, however, 
speeds up by the same amount, leading to an effective 
splitting in the rotating frame of 20. 
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IV. CONCLUSION 



We have shown that a time-dependent Lagrangian ap- 
proach provides physical insight into the nature of the 
low-energy modes of a vortex in a small trapped Bose- 
Einstein condensate. The anomalous mode is an in-phase 
circular motion of the vortex and the condensate with 
different amplitudes. In a rotating frame, the sense of 
this mode reverses from positive to negative when the 
external rotation frequency exceeds the metastable fre- 
quency jl4} . The dipole modes are rigid rotations of the 
vortex and the condensate together. In terms of the den- 
sity and phase fluctuations, all these results agree exactly 
with those from the Bogoliubov approach M . 
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